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ABSTRACT 


CaP” is analogous to the known interstellar anion CaN” with phosphorus replac¬ 
ing the nitrogen in a simple step down the periodic table. In this work, it is shown 
that CaP“ is likely to possess a dipole-bound excited state. It has been hypothesized 
and observationally supported that dipole-bound excited states are an avenue through 
which anions could be formed in the interstellar medium. Additionally, CaP” has a 
valence excited state that may lead to further stabilization of this molecule, and CaP” 
has a larger dipole moment than neutral CaP (~ 6 D vs. ^ 4 D). As such, CaP“ is 
probably a more detectable astromolecule than even its corresponding neutral radi¬ 
cal. Highly-accurate quantum chemical quartic force fields are also applied to CaP” 
and its singly substituted isotopologues in order to provide structures, vibrational 
frequencies, and spectroscopic constants that may aid in its detection. 


1 INTRODUCTION 


It is currently uncertain as to why anions are found 
in the interstellar medium (ISM) with six known to 
date: C 4 H-, CeH”, CsH”, CN”, CaN”, and CsN” 
(|Cemic haro_et_alJ|20Q2l; iBriinken et al.ll2007l: iRemiian et ahl 


2 OO 7 I: Thaddeu s et al. _ 20081: C ernicha ro et al.1 20081 : 

Agundez et al. I l2010l : ICordiner et al.l |2011^ . Since all of 


the known anions are closed-shell and, as a result, va¬ 
lence in nature, one of the leading theories about their 
existence co mes from the presen ce of a dipole-bound 
excited state (I Agundez et al.ll2008l ). If a neutral molecule 
possesses a large enough dipole moment (somewhere in 
excess of 2.0—2.5 D), it can bind an additional electron 
in a highly diffuse s-typ e orbital (iFermi fc Telleil 1 19471 : 
Coulson fc Walmslevl Il967|;_ Gntse^_&Adamowicz _ 1995blla : 


Compton et al.1119961: [Jordan fc Wang||2003l: ISimonsll200 


201l|). These loosely bound excited states dHammer et al.1 
2 OO 3 I ) can relax to a lower energy valence state if the 


molecule can support the additional electron in its valence 
electro n cloud. This phenomenon is known for so me small 
anions dLvkke et al.lll987l : lMullin et al.l[l9^ . Il993| ). and its 
reverse for electronic excitati on has even been hypothesized 
as a potential explanation (|Sarr 3 |2000l: fCordiner fc Sarrel 
I2OO7I : iFortenberry, Crawford fc Lee 2013l l for at least one 
of the diffuse interstellar bands, the series of visible to 
near-infrared absorption featur es largely un iformly observed 
towards many stellar objects dSarrd I2OO6I ) whose carriers 


have only very recently come into focus dCampbell et al.1 
I2OI5I) . As such, it is currently hypothesized that the 
large dipole moment of the neutral-radical captures an 
incident electron whereby the dipole-bound anion then 
decays down to the stable, valence ground ele ctronic state 
dAgiindez et al.l[2008l : IFortenberry et a^1l2013a^ . 

All of the known interstellar anions possess large enough 
dipole moments to undergo this process save for CN” 
which is likely formed by collisional or chemical pro cesses 
dbarsson. Geppert fc Nvmanll20l'3 : ICarelli et al.ll2014l ). The 
greater [CeH'J/ICsH] ratio as compared to the significantly 
smaller [C 4 H”]/[C 4 H] ratio is most eleg antly explained 
by th e dipole-bound formation hypothesis d Agundez et al .1 
I2OO8I) especially since C 4 H is more abundant than CeH, 
but Cell” is more abundant than C 4 H”. C 4 H must be ex¬ 
cited into its low-lying ^11 excited state in order to have a 
dipole moment large enough to capture the additional elec¬ 
tron. Cell is and strongly dipolar in its ground state 
meaning that capture of the dipole-bound electron is much 
more likely for the longer hydrocarbon chain. Additionally, 
C 2 II” has not been detec ted. Its corresponding neutral- 
radic al is very abundant dHeikkila. Johansson fc OlofssonI 
I1999I) . but the dipole moment of the unchanged species is 
not large enough to support a dipole-bound state. Finally, 
all of the corresponding neutral-radicals are known to ex- 
ist in the ISM for each of th e closed-shell anions detected 
llMcCarthv fc Thaddeusikooil ). 

However, no new anions have been detected in the 
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2 Fortenberry and Lukemire 


ISM since 2010 as reported in the literature llAgundez et alJ 
|2010|) . Several candidate anions have been studied quan¬ 
tum chemically in order to provide accurate predictions 
of their rovibrational nature. These include CH 2 CN~, c- 
and Z -CsH", and CCOH~ iFortenberrv, Crawford fc Lee 
2012 ; IFortenberrv et ^ I2013an2014al : IFortenberrv &: Lw 


20151 1. Another tantalizing potential interstellar anion 
is CsP”. It is analogous with the known CsN” an¬ 
ion, and several phosphorus-containing carbon chains 
CP, HCP, a nd CC P) have been detected in th e ISM 
Gu61in_^£^lJ_ 1990l: [A gudez. Cernicharo fc GuelinI l2007l : 


iHalfen, Clouthier fc Ziurvsll200^ . However, little is known 
about Z-CsP”. The linear neutral-radical is known to be the 
lowest energy isomer on the [C,C,C,P] potential energy sur- 
face and has a strongly dipo l ar (3. 7 97 D) ^11 ground stat e 
lldel Rio, Barrientos fc Largol 1 19961 : iMaclean et al.l l2008ll . 
Hence, this present work sets to provide both electronic as 
well as high-level rovibrational data on CaP” in order to 
assist in its potential detection. 

Quantum chemical studies of excited states of an¬ 
ions have recently shown that many anions may possess 
dipole-bound excited states, and s ome anion excited states 
may a ls o be valence in nature (IForten berrv &: Crawford 


2011blB iFortenb errvI 

20141 : IFortenberrv 


201F ; IFortenberrv. Morgan fc Envard 


2014 ; Theis et al.l 2015l l. Excited 


many species (ILvkke et al.l 

19871: iMullin et al.l 

1992, 

1993i: Grutter. Wvss & Maieil 

19991: iPachkov et al.l 

2001; 

Pino et al. 20021. 20041. but the number of dipole-bound 


excited state candidate anions has increased with the im¬ 
provement of computational tools. Additionally, high-level 
rovibrational computations have been able to produce 
rotational constants to within 15 MHz of experiment 
and vibrational frequencies to within 1.0 cm~^ of labo¬ 
ratory results in some cases ||Huan£j_T^^lor_&I^^2211 
Inostrozaj_Humi£_&_Jje^_l201l|_ Ftetenberr^^^_al 2012^ 


Hujmg^^jbrtenberry_^Led j20]jal : Zhao. Donev fc Linnarta 
20141 : Fortenberry et al. 2014bl l. Consequently, both ap¬ 
proaches are applied to CaP” in this work in order to 
elucidate more fully the chemical physics of this potentially 
astrochemically significant anion. 


2 COMPUTATIONAL DETAILS 


based on the optimized CCSD(T)/aug-cc-pVTZ geometry. 
Increasing the diffuseness of the basis set allows one to ob¬ 
serve the convergence of the excitation energy. Early con¬ 
vergence in the series indicates a valence excited state while 
late convergence is a sign that the excited state of interest 
is dipole-bound. This can them be corroborated by analy¬ 
sis of the orbitals involved from the most diffuse basis set. 

The double-zeta (n-aug-cc-pVDZ where n = 0 — 3) and 
triple zeta (m-aug-cc-pVTZ where m = 0 — 2) basis sets 
are employed. The vertical electron binding energy (eBE), 
the energy required to removed the electron from the sys¬ 
tem starting from the ground electronic state, is computed 
using the ionization pot ential form of EOM: EOMIP-CCSD 
dStanton fc Gau^ll994ll . 

The rovibrational co mputations ut iliz e the s cheme 
developed previous l y jHuang fc Led l2008l . l2009l : 

iHuang. Taylor fc Leel l2011^ that has met with much 
success as described above. The optimized CCSD(T)/aug- 
cc-pV5Z geometry is corrected for core-electron inclusion 
by adding the differences in the bond lengths between the 
CCSD(T) computations util izing the Martin-Taylor (MT) 
core correlating basis set llMartin fc Tavloil Il994^ with 
and without inclusion of the core electrons: Is for carbon 
and Is, 2s, and 2p for phosphorus. From this reference 
geometry, a quartic force field (QFF) is constructed. A 
QFF is fourth-order Taylor series expansi on of the potential 
portio n of the nuclear Hamiltonian (IFortenberrv et aP 
l2013bl l. It takes the form: 

U = - ^ Tij AiAj-b- ^ FijfeAiAj Afc-l-— ^ FijkiAiAjAkA.i, 

ij ijk ijkl 

( 1 ) 

where Ai are displacements and Fij... are the force con¬ 
stants. The displacements are defined from Ci as the termi¬ 


nal carbon atom opposite the phosphorus: 

5i(E+) = Ci-C2 (2) 

52(E+) = C2 - Cs (3) 

S3(E+) = C3-P (4) 

5'4(n,,) = /Cl - C 2 - C3 - y (5) 

5'5(n,,) = ZC2-C3-P-y (6) 

S6(nj,,) = zCi-C2-C3-x (7) 

S7(nj,,) = ZC2 - C3 - p - X. (8) 


All computations make use of coupled cluster the¬ 
ory at either the singles and doubles (CCSD) or sin- 
gles, doubles, and pertur b ative triples |CCSD(T)1 levels 
iRaghavacharl et al.l 198^ : ICrawford fc Schaefer Ilj I 2 OO 0 I : 
IShavitt fc Bartl^ 1200^ and the correll ation-consistent 
family of basis sets i n the c c-pVXZ from dPunnin iim. 
Peterson fc Dunnind 1 19951 : iKendall, Dunning fc HarrisonI 


1992h . The cc-pV('X-tdlZ basis sets are used for the phos¬ 


phorus atom, but this distinction should be understood by 
the reader to be implied for the remainder of the discussion. 
Additi onally, restricted H artree-Fock reference wavefunc- 
tions (IScheine r et al.ll l987^ and the PSI4 quantum chem¬ 
istry software | Turnev et al.ll20ia l are utilized except where 
noted. 

The electronically excited states are computed ver- 
tically using the equation of motion (EOM) formalism 
dStanton fc Bartl^ fl993l : iKrvlovI l2007h of CCSD and are 


A similar coordinate setup has been utilized pre- 
viously on the structurally-re lated Z-CsH^ cation 
([Huang, Fortenberry fc Leel l2013bl l. The step sizes for 
the bond length displacements are 0.005 A and 0.005 
radians for the bends. 

625 total displacements are needed to describe the 
CsP” QFF. At each displacement, seven quantum chemi¬ 
cal computations are undertaken at the CCSD(T) level. Ex¬ 
trapolation of the aug-cc-pVTZ, aug-cc-pVQZ, and aug-cc- 
pV5Z energies to th e complete basis set (CBS) limit via a 
three-point formula dMartln fc Leel[l996h provides the ref¬ 
erence energy for each point. Differences between the MT 
computed energies including and not including core elec¬ 
trons are added to the CBS energy. Furthermore, Douglas - 
Kroll scalar relativistic corrections dPouglas fc Krol]||l974fl 
availa ble in the MOLPRO 2010.1 quantum chemistry pro¬ 
gram dWerner et al.l [201^ are also added to the CBS en- 
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C'sP- 3 


ergy. The CcCR QFF (iFortenberrv et aDl201ll ~l is created 
from this approach where the CBS energy (“C”), core cor¬ 
relation (“cC”), and relativistic terms (“R”) provide for a 
very detailed description of the molecular potential energy 
surface. 

The energy points are fitted to the force constants via 
a least-squares formula. Refitting the surface provides the 
CcCR minimum equilibrium geometry as well as the proper 
force constants including zero value gradients. Transforma¬ 
tion of the force constants into Cartesian coordinates is 
accomplished via the INTDER l [Ajfen_fc cororkei^ l2005l l 
progra m. Rotational (iMillj 19^; IWatsonI 1977ll arid vibra¬ 
tional llPapousek &: Aliev! 119821 ) perturbation theory com¬ 
putations at second order (VPT2) produce spectroscopic 
constants and anharmonic vibrational frequencies, resp ec- 
tively, through the SPECTRO program llCaw et al.ll991| l. A 
2^'3 = vi type-1 Fermi resonance is included in SPECTRO’s 
VPT2 computations. Mpller-Plesset second-order electronic 
perturbat ion theory (MP2) computations with the 6-31-bG* 
basis set dMoller fc Plesseu 1 1934 iHehre. D itchfeld fc Pop lei 
I 1972 I ') in the Gaussian09 program ( Frisch et al.ll2009ll pro- 
vide the vibrational intensities. 


3 RESULTS AND DISCUSSION 

3.1 Vertical Electronically Excited States 

The ground electronic state C 3 P has a X term and 
a (core)7(T^8(j^9(j^l0(T^27r^ll (T^37r^ electronic configuration 
ddel Rfo, Barrientos fc Largol 1 1911^ 1. The anion’s additional 
electron occupies the 37 r orbital to complete the shell and 
give a valence X ground state term. The 37 r highest 
occupied molecular orbital (HOMO) is significant for the 
electronic properties of the anion. Computations of the first 
two vertically excited states of CsP” are given in Table [T] 
Both excitations are out of the 37 r HOMO. The 2 ^E+ state 
involves excitations into the Irr* orbital, while the 1 ^H state 
is produced by exciting the electron into a diffuse s-type or¬ 
bital. 

The 1 ^H state is a candidate for classification as 
a dipole-bound excited state. The excitation energy de¬ 
creases significantly as the spatial extent of the or¬ 
bitals increases, the excitations are, again, into high- 
diffuse orbitals as is necessary fo r a dipole-bound state 
IlSkurski. Gutowski fc Simons! I2OOOI : [Morgan fc Fortenberry! 
|2015^ . and the eBE and dipole-bound excitation energy are 
nearly coincident. Since the eBEs of di pole-bound states are 
known to be on the order of ~ 0.0 1 eV (iHammer et al.ir2003l : 
iDiken, Hammer fc Johnsonll2004 ). the latter point is a clear 
marker of dipole-bound behavior. Additionally, the dipole 
moment of the C 3 P X ^H state at the anion geometry is 
computed here by CCSD(T)/aug-cc-pVTZ to be 4.20 D, in 
line with previous computations llMaclean et al.l [20081) and 
large enough to support a dipole-bound state. 

The 2 ^E’^ state is actually valence. Its excitation is 
not affected by the increase in the diffuseness of the basis 
set and is well-below the eBE. It is the product of a tt —>■ tt* 
excitation which is similar to tt electron systems in polyun¬ 
saturated carbon chains, a commonly-used application of the 
particle-in-a-box quantum mechanical model system. Hence, 
CsP” and longer chains of this family will also likely possess 


a valence excited state. The CsP” 2 ^E"*" excitation energy 
is well converged to the 3.30 eV level for the double-zeta 
bases and similarly at 3.24 eV for the triple-zeta bases. 

Even though the intensities for the transitions are very 
small, GsP” is still an excellent candidate for dipole-bound, 
gas phase synthesis starting from the neutral-radical in the 
ISM or in circumstellar media. Capture of the electron by 
the neutral radical can proceed through the dipole-bound 
excited state producing the valence, closed-shell anion. The 
valence excited state provides an additional avenue for en¬ 
ergy dissipation. This may likely increase its interstellar life¬ 
time giving an exceptionally high [CsP^J/ICsP] ratio, po¬ 
tentially greater than that observed for [CeH'J/fCeH]. As a 
result, CsP” should be more readily observed than its cor¬ 
responding neutral radical since its electronic properties are 
more enhanced than the other known interstellar anions. In 
order to detect CsP”, information regarding the A^E'^ va¬ 
lence ground electronic state’s rotational, vibrational, and 
rovibrational behavior is needed and provided here for the 
first time from further quantum chemical computations. 


3.2 Rovibrational Analysis 

The fitting of the force constants needed to define the 
CcCR QFF and produce the subsequent spectroscopic data 
of ground state A^E"*" CsP” is very tight with a sum of 
residuals squared on the order of 10“^® a.u.^ The force con¬ 
stants for f-CsP” are given in Tableland follow the coordi¬ 
nate labeling given above. For instance, Fn corresponds to 
the second-derivative of the Ci—C 2 bond. Additionally, the 
harmonic, diagonal force constants are proportional to the 
bond strengths. Fn is the largest of these indicating that 
the Cl —C 2 bond is the strongest. This is corroborated by 
the Ra vibrationally averaged bond lengths given in Table 
[3] The ro(Ci—C 2 ) value is 1.259 055 A, the shortest bond 
length given indicating that this is the strongest bond. In¬ 
terestingly, the C 2 —C 3 bond strength from F 22 is roughly 
equal to the C 3 —P bond strength from F 33 . The phosphorus 
atom creates a longer bond length by default, but these force 
constants and the bond lengths show that this molecule is 
neither acetylenic nor cumulenic. It should be thought of 
as an intermediate between the two classifications with the 
Cl —C 2 bond more acetylenic and the rest of the molecule 
as cumulenic. 

In any case, the spectroscopic constants, dipole mo¬ 
ment, and vibrational frequencies are given in Table [3] along 
with the single ^®C substituted isotopologues. Phosphorus 
has only one stable isotope, ®^P. CsP” is strongly dipolar 
with a 6.22 D dipole moment as determined from center- 
of-mass computations. As a result, its rotational intensities 
should be brighter than the neutral-radical giving further 
evidence that the anion may be more easily observed in 
the ISM than the neutral. The rotational constants are not 
strongly affected by vibrational averaging with So equal to 
2798.0 MHz while is only 2.1 MHz greater at 2800.1 MHz. 
Substitution of the terminal Ci atom with ^®C reduces the 
rotational constants by 102.0 MHz. Substitution of the other 
carbon atoms is not as significant. The quartic (D) and sex- 
tic (H) distortion constants are exceptionally small for this 
molecule and may only show up in very high-resolution ex¬ 
periments. 

While ground-based telescopes such as the Atacama 
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4 Fortenberry and Lukemire 


Table 1. The C 3 P Vertical Excitation Energies (in eV), Ground State eBEs“ (in eV), and Oscillator Strengths.** 


Transition 

pVDZ 

apVDZ 

dapVDZ 

tapVDZ 

eBE 

/ 

pVTZ 

apVTZ 

dapVTZ 

2 is+ •«- 1 is+ 

3.52 

3.31 

3.30 

3.30 


2 X 10-"* 

3.35 

3.23 

3.24 

1 IR <- 1 is+ 

8.23 

5.04 

3.98 

3.81 

3.80 

8 X 10-3 

7.03 

4.89 

4.09 


“Computed with EOMIP-CCSD/t-aug-cc-pVDZ. 
^’The / values are for EOM-CCSD/d-aug-cc-pVTZ. 


Table 2. The CcCR C 3 P Force Constants Given in mdyn/A"'rad"*)“. 


Eli 

11.566 227 

F441 

-0.4428 

F2222 

219.70 

F5432 

-0.66 

Freii 

0.00 

F21 

0.596 174 

F442 

-0.5107 

F3111 

0.15 

F5433 

0.23 

F7e2i 

-1.13 

F22 

7.708 659 

F443 

0.0197 

F3211 

0.12 

F5444 

0.57 

F7S22 

0.77 

F31 

-0.380 686 

F541 

-0.0978 

F3221 

-0.38 

F5511 

-0.44 

F7S31 

-0.50 

F32 

0.689 835 

F542 

0.0236 

F3222 

0.66 

F5521 

0.65 

F7e32 

-0.19 

F33 

7.655 860 

F543 

-0.0500 

F3311 

-0.93 

F5522 

-0.88 

F7e33 

0.25 

F44 

0.209 473 

F551 

-0.0499 

F3321 

0.31 

F5531 

0.53 

F7S44 

-2.28 

F54 

-0.010 817 

F552 

-0.4968 

F3322 

0.11 

F5532 

1.80 

F7SS4 

3.47 

F55 

0.461 190 

F553 

-0.6287 

F3331 

1.77 

F5533 

0.62 

F7ess 

3.32 

Fee 

0.209 484 

Feel 

-0.4424 

F3332 

-0.64 

F5544 

0.71 

F7eee 

1.14 

Fre 

-0.010 826 

Fee2 

-0.5109 

F3333 

177.78 

F5554 

-0.03 

F7711 

-0.97 

F77 

0.461 209 

Fees 

0.0190 

F4411 

1.46 

F5555 

0.06 

F7721 

0.54 

Fill 

-71.2065 

F7ei 

-0.0952 

F4421 

-0.19 

Feeii 

1.09 

F7722 

-1.29 

F211 

-0.2476 

F7e2 

0.0233 

F4422 

0.61 

F6621 

-0.42 

F7731 

0.65 

F221 

-2.4219 

F7e3 

-0.0524 

F4431 

1.39 

F6622 

0.69 

F7732 

1.94 

F222 

-45.2600 

F771 

-0.0492 

F4432 

0.77 

Feesi 

1.53 

F7733 

-0.61 

F311 

-0.0455 

F772 

-0.4972 

F4433 

0.07 

F6632 

1.31 

F7744 

0.03 

F321 

0.7443 

F773 

-0.6289 

F4444 

-0.15 

Feess 

-0.38 

F77S4 

2.40 

F322 

-1.7436 

Fill! 

349.14 

F5411 

-0.05 

F6644 

1.89 

F77SS 

-0.68 

F331 

-0.0806 

F2111 

0.82 

F5421 

-0.66 

F6654 

-1.66 

F77es 

-0.08 

F332 

-0.8594 

F2211 

-3.24 

F5422 

-0.07 

Feess 

0.71 

F777e 

0.00 

F333 

-41.3790 

F2221 

4.66 

F5431 

0.64 

Feeee 

-0.25 

F7777 

-0.76 


“1 mdyn =10 ® N; The n and m values are exponents that correspond to the number of units from the type of modes present in the 

specific force constant. 


Large Millimeter/submillimeter (ALMA) array are excellent 
for observing molecular rotational spectra, the growth of 
space-based or airborne infrared telescopes such as the up¬ 
coming James Webb Space Telescope or the Stratospheric 
Observatory for Infrared Spectroscopy (SOFIA), respec¬ 
tively, has enhanced vibrational spectral analysis for astro¬ 
nomical applications. The vibrational frequencies of CsP” 
are given in the middle of Table [S] The descriptions of the 
fundamental vibrational modes are not as straightforward as 
one would assume for a Caov molecule. The first two modes 
are dehned here as the symmetric and antisymmetric C—C 
stretches for the three carbon atom portion of the molecule. 
The last two modes are the internal and external bends. The 
former mode is dominated by the C 2 and C 3 atoms moving 
opposite of one another relative to the central axis. The lat¬ 
ter bending vibrational motion is largely described by the 
Cl and P atoms moving in concert with one another rela¬ 
tive to the central axis. The carbon atom, naturally, covers 
a greater distance in this mode due to its smaller relative 
mass, but the two atoms are vibrating in concert. 

The CcCR QFF frequencies are not strongly anhar- 
monic. While anharmonic corrections are necessary for ac¬ 
curate descriptions of the vibrational modes, they do not 
shift as greatly from the harmonic values as is typical. 
Most of this behavior is probably due to the increased 
mass of the system largely from the phosphorus atom, but 


it is also from the tight bonding present in CsP”. An¬ 
ions are not often thought be tightly bound due to the 
increased electron-electron repulsion and likely increase in 
antibonding-character, but CsP” appears to be so bound 
from the data given. In fact, the C—C stretching frequen¬ 
cies in Cs P" are of the same order as thei r counterparts 
in C 3 H'*' liHuang. Fortenberry fc Led l2013bh . As a result, 
the vibrational frequencies provided here should be very 
good representations of physical reality. The frequencies ht 
SOFIA’S range very well, and the frequency at 1939.7 
cm“^ should be very bright since the double-harmonic MP2 
intensity is 483 km/mol. It is also worth mentioning that 
the MP2/6-31-I-G* harmonic frequency agree well with the 
CcCR harmonics . MP2 is know n to give Pauling points for 
some properties (ISherri Iil l2009lL but this correlation (< 5.1 
cm“^) is largely unprecedented but nicely coincidental for 
the significantly lower-level MP2 computations. 

In the age of ALMA, unidentihed lines (or “U-lines”) 
are becoming even more of an issue. Many of these are likely 
to be caused by rotational transitions of known molecules in 
vibrationally excited states. As a result, the vibrationally- 
averaged rotational constants of CsP” are also provided 
at the bottom of Table [3] for each fundamental mode. 
Vibrationally-excited modes associated with the bending 
frequencies increase the B values, while the stretching vibra¬ 
tional modes decrease them. substitution largely affects 


© 0000 RAS, MNRAS 000, 000-000 











C'sP- 5 


Table 3. The CsP” and Single Isotopologue CcCR Equilibrium and Zero-Point (Ra vibrationally-averaged) Geometries, Rotational 
Constants, Vibration-Rotation Interaction Constants, Quartic (D) and Sextic (H) Distortion Constants, CCSD(T)/aug-cc-pVTZ Dipole 
Moment, and Vibrational Frequencies and Intesnities“. 



CaP- 

i^cccp- 

ci^ccp- 

cci^cp- 

MP2/6-31+G* 

ro(Ci-C 2 ) (A) 

1.259 055 

1.259 200 

1.259 058 

1.259 065 


ro(C 2 -C 3 ) (A) 

1.345 935 

1.345 909 

1.345 920 

1.345 906 


ro(C 3 -P) (A) 

1.581 669 

1.581 612 

1.581 683 

1.581 652 


So (MHz) 

2798.0 

2696.0 

2769.6 

2798.0 


0 ® 1 (MHz) 

13.3 

13.2 

12.8 

12.9 


a® 2 (MHz) 

12.6 

11.8 

12.4 

12.4 


3 (MHz) 

4.3 

4.0 

4.2 

4.3 


a® 4 (MHz) 

-5.1 

-4.9 

-4.9 

-5.0 


a® 5 (MHz) 

-7.9 

-7.7 

-7.6 

-7.8 


re(Cl-C 2 ) (A) 

1.262 119 

- 

- 

- 


re(C 2 -C 3 ) (A) 

1.343 872 

- 

- 

- 


re(C 3 -P) (A) 

1.580 317 

- 

- 

- 


Be (MHz) 

2800.1 

2698.0 

2771.8 

2800.1 


De (kHz) 

0.202 

0.186 

0.199 

0.202 


He (/iHz) 

1.717 

1.654 

1.577 

1.717 


// ^ 

6.22 D 

- 

- 

- 


uji{cr) Symm. Carbon Stretch (cm“^) 

1979.2 

1959.0 

1930.9 

1971.8 

1984.3 (483) 

^ 2 ( 0 ') C 2 Shuttle (cm“^) 

1481.5 

1470.5 

1481.2 

1439.3 

1482.2 (93) 

a; 3 (cr) C 3 —P stretch (cm“^) 

692.0 

681.7 

686.6 

692.0 

689.2 (7) 

^ 4 ( 77 ) Internal Bend (cm~^) 

479.9 

479.4 

472.5 

469.6 

477.0 (4) 

a; 5 ( 7 r) External Bend (cm“^) 

194.8 

192.3 

192.4 

192.6 

195.5 (10) 

Harmonic Zero-Point (cm“^) 

2751.0 

2727.3 

2714.3 

2713.8 

2750.4 

Symm. Carbon Stretch (cm“^) 

1939.7 

1920.0 

1913.8 

1932.6 


^' 2 ( 0 ') C 2 Shuttle (cm“^) 

1456.5 

1446.0 

1456.0 

1415.5 


C 3 —P stretch (cm“^) 

684.0 

674.5 

678.2 

683.9 


z/ 4 ( 7 r) Internal Bend (cm“^) 

470.8 

470.4 

463.3 

461.3 


External Bend {cm~^) 

181.1 

178.8 

179.0 

179.4 


Zero-Point (cm~^) 

2730.0 

2706.7 

2698.4 

2693.8 


Bi (MHz) 

2784.7 

2682.8 

2756.8 

2785.1 


S 2 (MHz) 

2785.4 

2684.2 

2757.1 

2785.6 


S 3 (MHz) 

2793.7 

2692.0 

2765.3 

2793.8 


S 4 (MHz) 

2803.1 

2700.9 

2774.5 

2803.0 


S 5 (MHz) 

2805.9 

2703.7 

2777.2 

2805.8 



“Given in parentheses in km/mol beside the MP2/6-31-I-G* harmonic frequencies. 

^The coordinates used to generate the Born-Oppenheimer dipole moment component (in A with the centre-of-mass at the origin) are: 
Cl, 0.000000, 0.000000, -2.629136; C 2 , 0.000000, 0.000000, -1.367016; C 3 , 0.000000, 0.000000, -0.023143; P, 0.000000, 0.000000, 1.557174 


the rotational constants in similar ways between isotopo- 
logues as it does for Bq. Finally, a limited line list for the 
J ^ 9 rotational transitions of the v = Q and ^5 = 1 vi¬ 
brational states is provided in T able [4] in order assist with 
further analysis. Previous work l|Fortenberrv et ^ l2013al ) 
has shown that the rotational constants are probably within 
10 MHz of experiment putting the anticipated rotational 
constants within 2J x 10 MHz for each line predicted in 
Tabled and the relative intensities between the rotational 
transitions should be very close to experiment. 


4 CONCLUSIONS 

The electronic properties of CsP” are in line with those for 
the known interstellar anions indicating that it may be a 
viable molecule waiting to be detected in the ISM. Attach¬ 
ment of an additional electron to the corresponding neutral- 
radical via the anion’s dipole-bound excited state is a likely 


creation mechanism for this anion. Additionally, CaP” pos¬ 
sesses not only a valence ground electronic state but also 
a valence excited state. The second excited state should 
increase the molecule ’s lifetime. Since [C 6 H“]/[C 6 H] > 1 
llAgundez et ahlBoOSt) . it is likely that [C 3 P~]/[C 3 P] is as 
well since both radicals possess strongly dipolar ground elec¬ 
tronic states available for binding the extra electron leading 
to the observed valence anion. The valence excited state in 
CaP” should boost the [C 3 P“]/[C 3 P] ratio further due to 
the enhanced stability in the phosphorus anion. Further¬ 
more, the anion’s dipole moment is greater than that in the 
neutral-radical enhancing the anion’s rotational signal. As a 
result, CaP” may be a better candidate for detection than 
C 3 P. 

In order to assist in potential searches, the spectroscopic 
constants and vibrational frequencies for CsP” are provided 
in this work. Ground state CaP” is tightly bound, espe¬ 
cially considering that it is an anion, and not strongly anhar- 
monic. Hence, the CcCR approach utilized previously with 
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Table 4. The CsP” 150 K CcCR u = Q and 1/5 = 1 Rotational 
Limited Line List with AEj^i^j,J ^ 9 with Frequencies in 
MHz and Intensities in MHz per molecule. 


Mode 

j-fi 

J 

Frequency 

150K Intensity 

1 / = 0 

1 

0 

5596.0 

0.011 

u = 0 

2 

1 

11 192.0 

0.042 

u = 0 

3 

2 

16 788.0 

0.093 

u = 0 

4 

3 

22 384.0 

0.164 

u = 0 

5 

4 

27 980.0 

0.257 

u = 0 

6 

5 

33 575.9 

0.367 

u = 0 

7 

6 

39 171.8 

0.492 

u = 0 

8 

7 

44 767.7 

0.628 

u = 0 

9 

8 

50 363.5 

0.792 

u = 0 

10 

9 

55 959.3 

0.961 

1/5 = 1 

1 

0 

5611.8 

0.002 

1^5 = 1 

2 

1 

11 223.6 

0.006 

1/5 = 1 

3 

2 

16 835.4 

0.014 

1^5 = 1 

4 

3 

22 447.2 

0.025 

1^5 = 1 

5 

4 

28 058.9 

0.040 

t^5 = 1 

6 

5 

33 670.7 

0.057 

1^5 = 1 

7 

6 

39 282.4 

0.076 

1/5 = 1 

8 

7 

44 894.0 

0.097 

1/5 = 1 

9 

8 

50 505.7 

0.122 

1/5 = 1 

10 

9 

56 117.2 

0.148 


very good results should be similarly accurate in this sys¬ 
tem. The vi C—C stretch is very bright, and the vibration- 
rotation interaction constants are small. In snmmary, CaP” 
is an excellent interstellar molecular candidate, and the data 
provided in this work shonld assist in spectral searches for 
this anion whether in the laboratory or in the ISM. 
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